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In a recent work [Shao et al 2009 Phys. Rev. Lett. 108 018701], a nonconsensus opinion (NCO) model was 
proposed, where two opinions can stably coexist by forming clusters of agents holding the same opinion. The 
NCO model on lattices and several complex networks displays a phase transition behavior, which is character¬ 
ized by a large spanning cluster of nodes holding the same opinion appears when the initial fraction of nodes 
holding this opinion is above a certain critical value. In the NCO model, each agent will convert to its opposite 
opinion if there are more than half of agents holding the opposite opinion in its neighborhood. In this paper, we 
generalize the NCO model by assuming that each agent will change its opinion if the fraction of agents holding 
the opposite opinion in its neighborhood exceeds a threshold T (T > 0.5). We call this generalized model as the 
NCOT model. We apply the NCOT model on different network structures and study the formation of opinion 
clusters. We find that the NCOT model on lattices displays a continuous phase transition. For random graphs 
and scale-free networks, the NCOT model shows a discontinuous phase transition when the threshold is small 
and the average degree of the network is large, while in other cases the NCOT model displays a continuous 
phase transition. 
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I. INTRODUCTION 


The dynamics of opinion sharing and competing has become an active topic of recent research in statistical physics lH. 
One of the most successful methodologies used in opinion dynamics is agent-based modeling Q]. The idea is to construct the 
computational devices (known as agents with some properties) and then simulate them in parallel to model the real phenomena. 
In physics this technique can be traced back to Monte Carlo (MC) simulations ||2]. Beyond relevance as physics models, the 
ferromagnetic Ising model (l-llt], the XY model and the Potts model 10,01 can be seen as ^ent-based models for opinion 
dynamics. Other versions of opinion models have also been proposed, such as the Sznajd model 10, the majority rule model ifl^ 
the voter model id El, and the social impact model lUl . Some models display a disorder-order transition HHHtl, from 
a regime in which opinions are arbitrarily diverse to one in which most individuals hold the same opinion. Other models focus 
the emergence of a global consensus, in which all agents finally share the same opinion JEHU. 

It has been known that the formation of opinion clusters plays an important role in opinion dynamics An opinion 

cluster is defined as a connected component (subgraph) fully occupied by nodes holding the same opinion. Recently, Shao et 
al. proposed a nonconsensus opinion (NCO) model J^ in which each node adopts the majority opinion in its neighborhood at 
each time step. It was found that a large spanning cluster of nodes holding the same opinion appears when the initial fraction of 
nodes holding this opinion exceeds a certain threshold iddlldEll . Motivated by the NCO model, Li et al. proposed an inflexible 
contrarian opinion (ICO) model in which some agents never change their original opinion but may influence the opinions of 
others Ji^ . It was found that the threshold above which a large spanning cluster appears is increased with the fraction of 
inflexible contrarians. 

In both the NCO and ICO models, an agent will adopt the opinion that is held by more than half of neighbors. However, 
in many real-life situations, a quorum far larger than one half is necessary to pass a resolution. For example, a referendum to 
recall the president of the United States requires the support of two-thirds of the senators. Based on the above reasons, in this 
paper we generalize the NCO model by assuming that an agent will change its opinion when the fraction of agents holding the 
opposite opinion in its neighborhood exceeds a threshold T > 0.5. We call this generalized model as the NCOT model. When 
the threshold T = 0.5, the NCOT model recovers to the NCO model. When T = 1, the NCOT model becomes the standard 
percolation without opinion dynamics. Both the NCO and ICO models focus on the critical value for finite-size networks. By 
the standard finite-size scaling approach, we have obtained a critical point at which the phase transition takes place in the limit 
of infinite network size. It is interesting to And that, continuous or discontinuous phase transitions can arise in the NCOT model, 
depending on the value of T and the network structure. 
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FIG. 1: The normalized size of the largest cluster si and the second largest cluster S 2 as a function of / on a 1000 * 1000 square lattice. 
T = 0.7. Each curve is an average of 1000 different realizations. 


The paper is organized as follows. In Sec.|II] we introduce the NCOT model. In Sec. [Till we study the NCOT model on square 
lattices, random networks and scale-free networks, respectively. Finally, conclusions and discussions are presented in Sec. m 


II. A NONCONSENSUS OPINION MODEL WITH THE THRESHOLD (NCOT) 

In the NCOT model on networks, each node holds one of the binary opinions denoted by H-1 and —1. Initially, a fraction / 
of nodes with the opinion H-1 and 1 — / with the opinion —1 are selected at random. The neighborhood of node i is composed 
of node i and its nearest neighbors. At each time step, each node will convert to its opposite opinion, if the fraction of nodes 
holding the opposite opinion in its neighborhood exceeds a threshold T (T > 0.5). The system is considered to reach a stable 
state if no more changes occur. 


HI. MAIN RESULTS 


We focus on the formation of opinion clusters in the NCOT model. We denote by Si the size of the largest H-1 cluster and 
S 2 the size of the second largest H-1 cluster in the steady state. Then we define si = Si/N and S 2 = S 2 /N, where N is 
the network size. In the following, we carry out simulations systematically by employing the NCOT model on square lattices, 
random networks and scale-free networks respectively. 


A. The NCOT model on square lattices 

In this subsection, we study the NCOT model on an = L x L square lattice with periodic boundary conditions. Our 
extensive numerical simulations reveal that the phase transition of the NCOT model on square lattice can be roughly divided 
into three regimes (different universality classes of percolation) by tuning the parameter T, namely, T G [0.5, 0.6), [0.6,0.8) 
and [0.8,1]. In the same regime, the results are insensitive to the values of T. Particularly, it is noted that for 0.5 < T < 0.6, the 
corresponding phase transition pertains to the class of invasion percolation with trapping lldd]] : while for 0.8 < T < 1, it gives 
rise to a phase transition subjecting to regular site percolation. The percolation threshold is fc ~ 0.506 for T G [0.5,0.6) M 
and fc — 0.5927 for T G [0.8,1] respectively. For 0.6 < T < 0.8, the phase transition has some interesting features 
and has not been reported before, which will be a focus for this section. To be specific, in the following simulation, we choose 
T = 0.7. 

Figure[T]shows the normalized size of the largest cluster Si and the second largest cluster S 2 as a function of / when T = 0.7 
and N = 10®. We find that there exists a critical value fc{N), below which si approaches 0 and above which si continuously 
increases as / increases. At the critical value fc{N), S 2 displays a sharp peak, a characteristic of a second-order phase transi¬ 
tion Figure |2] shows the normalized size of the second largest cluster S 2 as a function of / for different values of N. From 






3 



/ 


FIG. 2: The normalized size of the second largest cluster S 2 as a function of / on a square lattice with different values of L. T = 0.7. Each 
curve is an average of 10000, 5000, 3000 and 1000 realizations for L — 100, 200, 400 and 1000, respectively. 



FIG. 3: (Color online) Log-log plot of (a) /c — fc{N), (b) the normalized size of the largest cluster si and (c) the susceptibility x, as a function 
of the system size N, respectively. The percolation threshold /c ~ 0.5492 for T = 0.7. In (a)-(c). The slopes of fitted lines are -0.36(1), 
-0.054(1) and 0.955(8) respectively. Each data point is an average of 10000, 5000, 4000, 3000, 2000 and 1000 realizations for L = 100, 200, 
300, 400, 600 and 1000, respectively. 


Fig.El we observe that the location of fc{N) changes with N. The percolation threshold fc{N) of a system of finite size N 
obeys the relation M 


fc{N) - ^ (1) 

where /c is the percolation threshold for a system of infinite size and iz is the correlation critical exponent. Then a simple 
linear fit (based on the maximization of the Pearson’s correlation coefficient) of fdN) vs. allows to simultaneously 

compute both values of /c and v At the percolation threshold fc, the normalized size of the largest cluster si and the 

susceptibility x = iV\/(si) — (si)^ vs. the system size N follow a power-law form: si ~ and x ~ 

From simulation results, we obtain ~ 0.5492 for T — 0.7. Figure [3 shows — fc{N), si and x as a function of N 
respectively. From Figs. |3a)-(c), we estimate the critical exponents l/zz = 0.36(1), fijv = 0.054(1) and 7 /iz = 0.955(8) as 
the best fit of the data points. 

For square lattice, our numerical simulations reveal that the critical point increases as T becomes larger. This is because 
for square lattice, the critical point fc is typically larger than 0.5. Thus in this region, a larger T hinders the transition of opinion 
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FIG. 4: The normalized size of the largest cluster si as a function of / for different values of T. The insets shows the normalized size of the 
second largest cluster S 2 versus / for different values of T. The ER network size N = 10® and the average degree (k) = 4. 



FIG. 5: The percolation threshold fc as a function of T on ER networks with the average degree (k) — 4. 


— 1 to +1 as it requires more +1 neighbors, reducing the final fraction of +1 nodes in the steady state, which in turn, will need 
a larger / for the spanning cluster of opinion +1 to emerge. 


B. The NCOT model on random networks 

In this subsection, we study the NCOT model on Erdos-Renyi (ER) random networks ER networks are characterized by 
a Poisson degree distribution with P{k) = {k)^/k\, where k is the degree of a node and (fc) is the average degree of the 

network. We perform simulations with different network sizes N. Each data point presented below is an average over 10000, 
8000, 6000, 4000, 3000, 2000 and 1000 different realizations for N = 10"^, 2 x 10"^, 5 x 10^, 10®, 2 x 10®, 5 x 10® and 10®, 
respectively. 

We have found that, when the average degree is small, the phase transition is continuous, while when the average degree is 
adequately large, depending on T, the transition can become discontinuous. 

To be specific, we have found that when the average degree (fc) = 4, the opinion percolation belongs to a continuous phase 
transition for the whole range of the parameter T. As shown in Eig. |4l for different values of T ranging from 0.5 to 1, the 
normalized size of the largest cluster si continuously increases with / and the normalized size of the second largest cluster S 2 
peaks at a certain critical value of /. Eigure |5] shows the percolation threshold fc as a function of T when the average degree 
(fc) = 4. The general trend is that as T increases, the critical threshold fc decreases. This can be understood that when T — 0.5, 
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FIG. 6: The normalized size of the largest cluster si as a function of / for (a) T = 0.5, (b) T = 0.6 and (c) T = 0.8. The ER network size 
N = 10® and the average degree (fe) = 10. 


a node switches its opinion when half of its neighbors (including itself) has the opposite opinion, therefore nodes tend to form 
clusters with the same opinions. But since T is small, the cluster is not compact, but rather sparsely connected. In the region 
where the initial fraction of nodes with +1 opinion / < 0.5, nodes with -1 are majority, the node is more likely to be surrounded 
by neighbors with opposite opinions. Thus if T is small, a node with +1 is more likely to switch to -1, therefore the effective 
fraction of nodes with +1 becomes smaller, making it more difficult for a spanning cluster to emerge. As T increases, it is harder 
for a node to switch its opinion (requires more neighbors with opposite opinions), leading to a larger effective fraction of nodes 
with +1, and therefore a higher probability for a spanning cluster to emerge. Note that the effective fraction of nodes with +1 
will still be smaller than / since nodes with -1 are still majorities therefore more nodes with +1 will be switched to -1 than the 
opposite process. 

A surprising phenomenon is that there exists an optimal value of T (about 0.7) leading to the minimum of fc- This could 
be a result of higher order interaction between the network topology and the opinion dynamics. For example, although a larger 
T means more difficult for a node to switch its opinion, but once the condition is satisfied and this node switches, the cluster 
grows and it will be more compact, and it will be more resistive to changes caused by outside nodes, leading to a nonmonotonic 
behavior of fc versus T. The nonmonotonic relation between fc and T can also confirmed in Fig.lH From Fig.lH we can see 
that the critical value of / that corresponds to the peak of S 2 is the smallest when T = 0.7. 

When the average degree (fc) is adequately large (e.g., (fc) = 10), the opinion percolation can display a continuous or a 
discontinuous phase transition, depending on the value of T. Figure | 6 ] shows the normalized size of the largest cluster si as a 
function of / for different values of T when (fc) = 10. From Figs.lhja) and (b), we observe that there exist two abrupt transition 
points when T is small (e.g., T — 0.5 and T = 0.6). At the first abrupt transition point denoted by fc, si jumps from zero 
to a finite value. At the second abrupt transition point denoted by /*, si jumps from a finite value to one. Between the two 
abrupt transition points, Si continuously increases with /. However, for the large value of T (e.g., T = 0.8), si approaches 
zero continuously as / is decreased from 1 to 0 [see Fig.| 6 lc)]. These numerical results indicate that the behavior of si versus 
/ could be a discontinuous phase transition for the small value of T and a continuous phase transition for the large value of T. 
From Fig.| 6 ] it is seen that the system has a symmetry, i.e., the curves are unchanged if / 1 — / and Si 1 — Si. This is 

because, in this case, the nodes are densely connected, therefore one can neglect the small clusters for either +1 or —1 opinions. 
Denote sj as the normalized size for the largest cluster for —1, we have si + s* ~ 1. Since the dynamics for the evolution of 
+1 into —1 and vice versa are the same, this imposes a duality between +1 and —1 states. Therefore, for a given / after the 
opinion dynamics reach the steady state the normalized size of the largest cluster for +1 is si, is actually the same process for 
— 1 opinions with a given initial fraction 1 — / and a largest cluster with normalized size sj(/) ft! 1 — si(/). But since the 
dynamics for +1 and —1 are the same, when the initial fraction for +1 is 1 — /, the largest cluster for +1 will also be s*(/), 
which is approximately 1 — Si(/). 

Figure |7] shows that si as a function of / for different values of the network size N when (fc) = 10 and T = 0.6. One can 
see that all the curves intersect at one point (the phase transition point). In Fig.[^ we investigate the number of iterations (NOI), 
which is the number of time steps needed to reach the steady state, as a function of / when (fc) = 10 and T = 0.6. Note that 
NOI characterizes the long range correlation, i.e., if the correlation is local, the system will quickly settle down to the steady 
state; while if there exist long range correlations, it needs more iterations to reach the steady state since each status change for a 
node has wider impacts. We can observe that the NOI exhibits two symmetric peaks. According to Ref. ll4^ . in a discontinuous 
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FIG. 7: The normalized size of the largest cluster si as a function of / for different values of the network size N. The average degree of ER 
networks is (k) = 10 and T = 0.6. The initial fraction of nodes with the +1 opinion / is around fc for (a) and / is around for (b). 



/ 


FIG. 8: The number of iterations (NOI) as a function of / on ER networks. The network size N = 10®, the average degree (fe) = 10 and 

r = o.6. 


phase transition, the location of the peak of the NOI determines the critical threshold of the transition. In Fig. 0 the location 
of the left peak determines the critical threshold fc below which si = 0 and the right peak determines the critical threshold /* 
above which si = 1. From simulation results shown in Figs. 7 and 8 , for (k) = 10 and T = 0.6, we obtain fc — 0.4408 and 
/* ~ 0.5592. It is noted that /c + /* = 1, which is consistent with the previous symmetry analysis. 

To further classify the transition class, we carry out finite size scaling analysis. Figure |9] shows the normalized size of the 
largest cluster si, the susceptibility x and the number of iterations (NOI), as a function of the system size N at the discontinuous 
transition points fc and /*, respectively. Figure|9ta) shows that si scales as with 13/v Pi Q, indicating a discontinuous 

phase transition. Figurej^b) illustrates that x scales as with 7 / 1 / Ri 1. Figure|9jc) shows that NOI scales as with 

5 ~ 0.28, consistent with the theoretical result 1/4 

Figure [To] shows the percolation threshold fc as a function of T when the average degree {k) = 10. One can see that fc 
decreases to 0.1 as T increases. There exists a certain critical value Tc (between 0.7 and 0.75), below which the phase transition 
is discontinuous while above which the phase transition becomes continuous. The development to the first-order-like phase 
transition is originated to the consolidation of clusters with the same opinions Two key factors are needed. The first 

one is the network topology, which needs to be dense enough in order to promote clustering process during the evolution of 
the opinions S. The second one is the opinion dynamics, where T should be small to lower the barrier for a node to switch 
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FIG. 9: Log-log plot of (a) the normalized size of the largest cluster si, (b) the susceptibility x and (c) the number of iterations (NOI), as a 
function of the system size N, respectively. In (a)-(c). The slopes of fitted lines are about 0, 1 and 0.28 respectively. The average degree of ER 
networks (k) = 10 and T = 0.6. The phase transition points are /c ~ 0.4408 and /* ~ 0.5592. 



FIG. 10: The percolation threshold fc as a function of T on ER networks with the average degree (k) = 10. The phase transition is 
discontinuous for small values of T (filled circles), while it becomes continuous for large values of T (empty circles). 


its opinion to facilitate the formation of clusters. Therefore, in the region when (k) is large and T is small, one could expect 
disrupt emerging of spanning clusters; while in the opposite case when (k) is small and T is large, one would expect continuous 
transitions. 


C. The NCOT model on scale-free networks 

In this subsection, we study the NCOT model on Barabasi-Albert scale-free networks (BA) if^ . BA networks are character¬ 
ized by a power-law degree distribution with P{k) ^ k~^. We perform simulations with different network sizes N. Each data 
point is an average over 10000, 6000, 4000, 3000, 2000 and 1000 different realizations for N = 10^, 2 x 10^, 5 x 10"*, 10®, 
2 X 10® and 5 x 10®, respectively. 

Figure [TT] shows the normalized size of the largest cluster si as a function of / for different values of T. We can see that, for 
(fc) = 4 (Fig. fTTl a)). Si approaches zero continuously as / decreases for different values of T, indicating a continuous phase 
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FIG. 11: The normalized size of the largest cluster si as a function of / for different values of T. The average degree of BA networks is (a) 
(fe) = 4 and (h) (fc) = 12, respectively. The network size N = 5 x 10®. 



FIG. 12: The percolation threshold /c as a function of T for different values of the average degree (k) of BA networks. Filled circles denote 
that the phase transition is discontinuous, while empty circles and empty triangles represent that the transition is continuous. 


transition. While for (k) = 12 (Fig. fTTl b)). there exists two abrupt transition points when T is small (e.g., T = 0.6), and the 
transition becomes continuous for a larger T {T = 0.8). At the hrst abrupt transition point denoted by fc, si jumps from zero 
to a finite value. At the second abrupt transition point denoted by /*, si jumps from a bnite value to one. We have checked that 
fc + fc = which is the same as that in ER networks. 

Figure [12] shows the percolation threshold fc as a function of T for different values of (k). One can see that fc decreases as 
the increase of T. When the average degree (fc) is small (e.g., (fc) = 4), the percolation belongs to a continuous phase transition 
for all the values of T. When (fc) is large (e.g., (fc) = 12), the percolation behaves a discontinuous transition for the small values 
of T while it displays a continuous phase transition when T is large. 


IV. CONCLUSIONS AND DISCUSSIONS 

In conclusion, we have proposed a generalized nonconsensus opinion model in which an agent changes its opinion when the 
fraction of nodes holding the opposite opinion in its neighborhood exceeds a threshold T {T > 0.5). We apply the model on 
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various network structures to study the formation of opinion clusters. It is found that the behavior of the normalized size of 
the largest cluster versus the initial concentration of nodes holding the same opinion can display a discontinuous or continuous 
phase transition. For regular lattices, the phase transition is continuous independent of T. For complex networks such as random 
networks and scale-free networks, if the average degree is small, then the phase transition is continuous, regardless of the value 
of T. For complex networks with the large average degree, the phase transition is continuous when T is large but it becomes 
discontinuous when T is small. Particulary, there exists two symmetric critical values in the case of the discontinuous phase 
transition. The studied opinion disappears below the first critical value while it takes over the whole population above the second 
critical value. 

We also study the relationship between T and the percolation threshold /c above which a large spanning cluster of nodes 
holding the same opinion appears. Note that when T = 1, the phase transition of our model reverts to the regular site percolation 
and the opinion percolation threshold is equal to that of site percolation. For square lattices, the percolation threshold increases 
as T increases from 0.5 to 1. For Erdos-Renyi random networks with the small values of the average degree, there exists an 
optimal value of T, leading to the minimum fc- For Erdos-Renyi random networks with the large values of the average degree 
or Barabasi-Albert scale-free networks, fc decreases as the increase of T. 
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